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Abstract 

In this paper, we consider a scenario where a source node wishes to broadcast two confidential messages for two respective 
receivers via a Gaussian MIMO broadcast channel. A wire-tapper also receives the transmitted signal via another MIMO channel. 
First we assumed that the channels are degraded and the wire-tapper has the worst channel. We establish the capacity region 
of this scenario. Our achievability scheme is a combination of the superposition of Gaussian codes and randomization within 
the layers which we will refer to as Secret Superposition Coding. For the outerbound, we use the notion of enhanced channel 
to show that the secret superposition of Gaussian codes is optimal. We show that we only need to enhance the channels of the 
legitimate receivers, and the channel of the eavesdropper remains unchanged. Then we extend the result of the degraded case to 
non-degraded case. We show that the secret superposition of Gaussian codes along with successive decoding cannot work when 
the channels are not degraded, we develop a Secret Dirty Paper Coding (SDPC) scheme and show that SDPC is optimal for this 
channel. Finally, we investigate practical characterizations for the specific scenario in which the transmitter and the eavesdropper 
have multiple antennas, while both intended receivers have a single antenna. We characterize the secrecy capacity region in terms 
of generalized eigenvalues of the receivers channel and the eavesdropper channel. We refer to this configuration as the MISOME 
case. In high SNR we show that the capacity region is a convex closure of two rectangular regions. 



I. Introduction 

Recently there has been significant research conducted in both theoretical and practical aspects of wireless communication 
systems with Multiple-Input Multiple-Output (MIMO) antennas. Most works have focused on the role of MIMO in enhancing 
the throughput and robustness. In this work, however, we focus on the role of such multiple antennas in enhancing wireless 
security. 

The information-theoretic single user secure communication problem was first characterized by Wyner in [1]. Wyner 
considered a scenario in which a wire-tapper receives the transmitted signal over a degraded channel with respect to the 
legitimate receiver's channel. He measured the level of ignorance at the eavesdropper by its equivocation and characterized the 
capacity-equivocation region. Wyner's work was then extended to the general broadcast channel with confidential messages 
by Csiszar and Korner [2]. They considered transmitting confidential information to the legitimate receiver while transmitting 
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a common information to both the legitimate receiver and the wire-tapper. They established a capacity-equivocation region of 
this channel. The secrecy capacity for the Gaussian wire-tap channel was characterized by Leung- Yan-Cheong in [3]. 

The Gaussian MIMO wire-tap channel has recently been considered by Khisti and Wornell in [4], [5]. Finding the optimal 
distribution, which maximizes the secrecy capacity for this channel is a nonconvex problem. Khisti and Wornell, however, 
followed an indirect approach to evaluate the secrecy capacity of Csiszar and Korner. They used a genie-aided upper bound 
and characterized the secrecy capacity as the saddle-value of a min-max problem to show that Gaussian distribution is optimal. 
Motivated by the broadcast nature of the wireless communication systems, we considered the secure broadcast channel with an 
external eavesdropper in [6], [7] and characterized the secrecy capacity region of the degraded broadcast channel and showed 
that the secret superposition coding is optimal. Parallel and independent with our work of [6], [7], Ekrem et. al. in [8], [9] 
established the secrecy capacity region of the degraded broadcast channel with an external eavesdropper. The problem of 
Gaussian MIMO broadcast channel without an external eavesdropper is also solved by Lui. et. al. in [10]— [12]. 

The capacity region of the conventional Gaussian MIMO broadcast channel is studied in [13] by Weingarten et al. The notion 
of an enhanced broadcast channel is introduced in this work and is used jointly with entropy power inequality to characterize 
the capacity region of the degraded vector Gaussian broadcast channel. They showed that the superposition of Gaussian codes 
is optimal for the degraded vector Gaussian broadcast channel and that dirty-paper coding is optimal for the nondegraded case. 

In the conference version of this paper (see [14]), we established the secrecy capacity region of the degraded vector Gaussian 
broadcast channel. Our achievability scheme, was a combination of the superposition of Gaussian codes and randomization 
within the layers which we refereed to as Secret Superposition Coding. For the outerbound, we used the notion of enhanced 
channel to show that the secret superposition of Gaussian codes is optimal. In this paper, we aim to characterize the secrecy 
capacity region of a general secure Gaussian MIMO broadcast channel. Our achievability scheme is a combination of the dirty 
paper coding of Gaussian codes and randomization within the layers. To prove the converse, we use the notion of enhanced 
channel and show that the secret dirty paper coding of Gaussian codes is optimal. We investigate practical characterizations 
for the specific scenario in which the transmitter and the eavesdropper have multiple antennas, while both intended receivers 
have a single antenna. This model is motivated when a base station wishes to broadcast secure information for small mobile 
units. In this scenario small mobile units have single antenna while the base station and the eavesdropper can afford multiple 
antennas. We characterize the secrecy capacity region in terms of generalized eigenvalues of the receivers channel and the 
eavesdropper channel. We refer to this configuration as the MISOME case. In high SNR we show that the capacity region is a 
convex closure of two rectangular regions. 

Parallel with our work, Ekrem et. al [15] and Liu et. al. [16], [17], independently considered the secure MIMO broadcast 
channel and established its capacity region. Ekrem et. al. used the relationships between the minimum-mean-square-error 
and the mutual information, and equivalently, the relationships between the Fisher information and the differential entropy to 
provide the converse proof. Liu et. al. considered the vector Gaussian MIMO broadcast channel with and without an external 
eavesdropper. They presented a vector generalization of Costa's Entropy Power Inequality to provide their converse proof. In 
our proof, however, we enhance the channels properly and show that the enhanced channels are proportional. We then use 
the proportionality characteristic to provide the converse proof. The rest of the paper is organized as follows. In section II we 
introduce some preliminaries. In section III, we establish the secrecy capacity region of the degraded vector Gaussian broadcast 
channel. We extend our results to non-degraded and non vector case in section IV. In Section V, we investigate the MISOME 
case. Section VI concludes the paper. 
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Fig. 1. Secure Gaussian MIMO Broadcast Channel 
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II. Preliminaries 

Consider a Secure Gaussian Multiple-Input Multiple-Output Broadcast Channel (SGMBC) as depicted in Fig. Q] In this 
confidential setting, the transmitter wishes to send two independent messages (Wi, Wb) to the respective receivers in n uses 
of the channel and prevent the eavesdropper from having any information about the messages. At a specific time, the signals 
received by the destinations and the eavesdropper are given by 



yi = Hix + ni, 
y 2 = H 2 x + n 2 , 

Z = H3X + II3, 



(1) 



where 

• x is a real input vector of size t X 1 under an input co variance constraint. We require that E[xx T ] < S for a positive 
semi-definite matrix S X 0. Here,^;, ^. >-, and X represent partial ordering between symmetric matrices where B X A 
means that (B — A) is a positive semi-definite matrix. 

• yi> Y2, and z are real output vectors which are received by the destinations and the eavesdropper respectively. These are 
vectors of size n x 1, ri x 1, and ra x 1, respectively. 

• Hi, H2, and H3 are fixed, real gain matrices which model the channel gains between the transmitter and the receivers. 
These are matrices of size 7*1 x t, r 2 x t, and x t respectively. The channel state information is assumed to be known 
perfectly at the transmitter and at all receivers. 

• ni, n 2 and 113 are real Gaussian random vectors with zero means and covariance matrices Ni = ^[nin^] >- 0, 
N 2 = -B[n 2 ii2 T ] X 0, and N 3 = £[n 3 ii3 T ] X respectively. 

Let W\ and W2 denote the the message indices of user 1 and user 2, respectively. Furthermore, let X, Y\, Y2, and Z denote 
the random channel input and random channel outputs matrices over a block of n samples. Let V\, V2, and V3 denote the 
additive noises of the channels. Thus, 



Yi=H!X + V u 

Y 2 = h 2 T + V 2 , 
1 = H3X + F3. 



(2) 



Note that Vi is an x n random matrix and Hj is an r, x t deterministic matrix where i — 1, 2, 3. The columns of Vi are 
independent Gaussian random vectors with covariance matrices Ni for i — 1,2,3. In addition Vi is independent of X, W% 
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and W%. A ((2 nRl ,2 nR2 ),n) code for the above channel consists of a stochastic encoder 

/ : ({1, 2, ...,2 nRl ] x {1, 2, 2 nR *}) -» (3) 

and two decoders, 

gi :%^{l,2,...,2 nR i}, (4) 

and 

5 2 :^2^{l,2,...,2" fl2 }- (5) 

where a script letter with double overline denotes the finite alphabet of a random vector. The average probability of error is 
defined as the probability that the decoded messages are not equal to the transmitted messages; that is, 

Pj n ) = P( gi (Fi ) ^ Wi U 52 (Fa ) # W 2 ) . (6) 

The secrecy levels of confidential messages W\ and Wa are measured at the eavesdropper in terms of equivocation rates, 
which are defined as follows. 

Definition 1: The equivocation rates R e i, R e2 and R e \ 2 for the secure broadcast channel are: 

R el = -H(W 1 ffl, (7) 
n 

R e2 = ~H{W 2 \Z), 

n 

R eU = -H{W U W 2 \2). 

n 

The perfect secrecy rates R\ and R 2 are the amount of information that can be sent to the legitimate receivers both reliably 
and confidentially. 

Definition 2: A secrecy rate pair (Ri,R 2 ) is said to be achievable if for any e > 0, ei > 0, e 2 > 0,e3 > 0, there exists a 
sequence of ((2 nRl , 2 nR2 ) 1 n) codes, such that for sufficiently large n, 

P( n) <e, (8) 
R e i>Ri-ei, (9) 
R e2 >R 2 ~e 2 , (10) 

R e i 2 >Ri + R 2 -e 3 . (11) 
In the above definition, the first condition concerns the reliability, while the other conditions guarantee perfect secrecy for each 
individual message and both messages as well. The model presented in (Q~|) is SGMBC. However, we will initially consider 
two subclasses of this channel and then generalize our results for the SGMBC. 

The first subclass that we will consider is the Secure Aligned Degraded MIMO Broadcast Channel (SADBC). The MIMO 
broadcast channel of (Q} is said to be aligned if the number of transmit antennas is equal to the number of receive antennas at 
each of the users and the eavesdropper (t = r\ = r 2 = rs) and the gain matrices are all identity matrices (Hi = H2 = H3 = I). 
Furthermore, if the additive noise vectors' covariance matrices are ordered such that -< Ni ^ N2 d N3, then the channel 
is SADBC. 

The second subclass we consider is a generalization of the SADBC. The MIMO broadcast channel of (HJ is said to be Secure 
Aligned MIMO Broadcast Channel (SAMBC) if it is aligned and not necessarily degraded. In other words, the additive noise 
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vector covariance matrices are not necessarily ordered. A time sample of an SAMBC is given by the following expressions, 

yi=x + m, (12) 
y 2 = x + n 2 , 
z = x + n 3 , 

where, yi, y 2 , z, x are real vectors of size t x 1 and nx, n 2 , and n3 are independent and real Gaussian noise vectors such 
that Ni = E[niiii T ] >- txt for i = 1, 2, 3. 

III. The Capacity Region of The SADBC 

In this section, we characterize the capacity region of the SADBC. In [6], we considered the degraded broadcast channel 
with confidential messages and establish its secrecy capacity region. 

Theorem 1: The capacity region for transmitting independent secret messages over the degraded broadcast channel is the 
convex hull of the closure of all (i?i,i?2) satisfying 

Ri<I{X\Yx\U)-I{X\Z\\J\ (13) 
R 2 <I(U;Y 2 )-I(U;Z). (14) 

for some joint distribution P(u)P(x\u)P(yi, y 2 , z\x). 

Proof: Our achievable coding scheme is based on Cover's superposition scheme and random binning. We refer to this 
scheme as the Secret Superposition Scheme. In this scheme, randomization in the first layer increases the secrecy rate of the 
second layer. Our converse proof is based on a combination of the converse proof of the conventional degraded broadcast 
channel and Csiszar Lemma. Please see [6], [7] for details. ■ 
Note that evaluating ( fT3T > and ( fl4l i involves solving a functional, nonconvex optimization problem. Usually nontrivial techniques 
and strong inequalities are used to solve optimization problems of this type. Indeed, for the single antenna case, [18], [19] 
successfully evaluated the capacity expression of ([T3l and (TBI) . Liu et al. in [20] evaluated the capacity expression of MIMO 
wire-tap channel by using the channel enhancement method. In the following section, we state and prove our result for the 
capacity region of SADBC. 

First, we define the achievable rate region due to Gaussian codebook under a covariance matrix constraint S >z 0. The 
achievability scheme of Theorem [T] is the secret superposition of Gaussian codes and successive decoding at the first receiver. 
According to the above theorem, for any covariance matrix input constraint S and two semi-definite matrices Bi >; and 
B 2 y such that Bi + B 2 < S, it is possible to achieve the following rates, 



i?f(B 1 , 2 ,N 1 , 2 , 3 ) = i 
fl2 G (Bx, 2 ,N 1)2 , 3 ) = ^ 



log I Ni 1 (Bx + Nx) | - ~ log | N 3 1 (Bx + N 3 ) | 
'. (Bi+Ba+Nsl 1 IBx+Ba + Nal 

lo § — F5 — n^n o log ■ 



|Bx+N 2 | 2 ° iBi+Nal 
The Gaussian rate region of SADBC is defined as follows. 

Definition 3: Let S be a positive semi-definite matrix. Then, the Gaussian rate region of SADBC under a covariance matrix 
constraint S is given by 

' (i^(B 1 , a ,N 1>a) 3),i^(B 1 , a ,N 1>a)3 )) | 
s.t S - (Bi + B 2 ) t0, B k h 0, k = 1,2 



ft G (S,N l!2 , 3 ) - 
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We will show that 1Z G (S, Ni i2i3 ) is the capacity region of the SADBC. Before that, certain preliminaries need to be addressed. 
We begin by characterizing the boundary of the Gaussian rate region. 

Remark 1: Note that in characterizing the capacity region of the conventional Gaussian MIMO broadcast channel Weingarten 
et al. [13] proved that on the boundary of the above region we have Bi + B2 = S which maximizes the rate R 2 . In our 
argument, however, the boundary is not characterized with this equality as rate R 2 may decreases by increasing Bi + B2. 

Definition 4: The rate vector R* = (R\, R 2 ) is said to be an optimal Gaussian rate vector under the co variance matrix S, if 
R* e ft G (S,Ni j2 ,3) and if there is no other rate vector R'* = (i?i,i? 2 ) € fc G (S,Ni j2i3 ) such that R[ > R 1 and R' 2 > R 2 
where at least one of the inequalities is strict. The set of positive semi-definite matrices (BJ,B 2 ) such that BJ + B 2 < S is 
said to be realizing matrices of an optimal Gaussian rate vector if the rate vector (i? G (B£ 2 , Ni i2i3 ), i? G (BJ 2 , Ni i2 ,3)) is 
an optimal Gaussian rate vector. 

In general, there is no known closed form solution for the realizing matrices of an optimal Gaussian rate vector. Note that 
finding an optimal Gaussian rate vector once again, involves solving a nonconvex optimization problem. The realizing matrices 
of an optimal Gaussian rate vector, BJ, B 2 are the solution of the following optimization problem: 

max fl G (B li2 ,N 1;2;3 )+/^ G (B 1;2 ,N li2i3 ) (15) 

(Bi,B 2 ) 

s.tBi^O, B 2 ^0, Bi+B 2 ^S, 

where fi > 1. Next, we define a class of enhanced channel. The enhanced channel has some fundamental properties which 
help us to characterize the secrecy capacity region. We will discuss its properties later on. 

Definition 5: A SADBC with noise covariance matrices (N 1 ,N 2 ,N 3 ) is an enhanced version of another SADBC with 
noise covariance matrices (Ni,N 2 ,N 3 ) if 

N'^Nx, N 2 ^N 2 , N 3 = N 3 , N;^N 2 . (16) 
Obviously, the capacity region of the enhanced version contains the capacity region of the original channel. Note that in 
characterizing the capacity region of the conventional Gaussian MIMO broadcast channel, all channels must be enhanced by 
reducing the noise covariance matrices. In our scheme, however, we only enhance the channels for the legitimate receivers and 
the channel of the eavesdropper remains unchanged. This is due to the fact that the capacity region of the enhanced channel 
must contain the original capacity region. Reducing the noise covariance matrix of the eavesdropper's channel, however, may 
reduce the secrecy capacity region. The following theorem connects the definitions of the optimal Gaussian rate vector and 
the enhanced channel. 

Theorem 2: Consider a SADBC with positive definite noise covariance matrices (Ni, N 2 , N 3 ). Let B^ and B 2 be realizing 
matrices of an optimal Gaussian rate vector under a transmit covariance matrix constraint S >- 0. There then exists an enhanced 
SADBC with noise covariance matrices (N 1 ,N 2 ,N 3 ) that the following properties hold. 

1) Enhancement: 

n;^n 15 n 2 ^n 2 , n 3 = n 3 , n;^n 2 , 

2) Proportionality: 

There exists an a > and a matrix A such that 
(I-A)(B* 1 +N' 1 )=aA(Bi+N 3 ), 

3) Rate and optimality preservation: 

^f(Bi,2! Ni i2)3 ) = -R G (Bi ;2 , N 12 3 ) V/c = 1,2, furthermore, B£ and B 2 are realizing matrices of an optimal 
Gaussian rate vector in the enhanced channel. 
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Proof: The realizing matrices B\ and B 2 are the solution of the optimization problem of ( TT3T ). Using Lagrange Multiplier 
method, this constraint optimization problem is equivalent to the following unconditional optimization problem: 

max Rf (Bi, a , N 1)2 , 3 ) + fiR§(B ly2 , N ll2)3 ) + Tr{BiOi} 

(Bi ,B 2 ) 

+Tr{B 2 2 } + Tr{(S — Bi — B 2 )0 3 }, 

where Oi, 2 , and O3 are positive semi-definite t x t matrices such that TrjBJOi} = 0, Tr{B 2 2 } — 0, and 
Tr{(S - BJ - B 2 )0 3 } = 0. As all B£, k = 1,2, O i; i = 1,2,3, and S B£ B 2 are positive semi-definite matrices, 
then we must have B£Ok = 0, k = 1, 2 and (S — B£ — B 2 )0 3 = 0. According to the necessary KKT conditions, and after 
some manipulations we have: 

(bi + Ni)- 1 + Oi - i)(b; + N3)- 1 + Ox = m (b; + n,)- 1 + o 2 , (17) 

M( B i + B 2 + N,)- 1 + 2 = M (BJ + B 2 + Na)- 1 + O s . (18) 
We choose the noise covariance matrices of the enhanced SADBC as the following: 

N' 1= (Nr'+Oi)" 1 , (19) 
N 2 = ((B\ + N?)- 1 + ~0 2 ) - B* l7 
N 3 = N 3 . 

As Oi ^ and 2 y 0, then the above choice has the enhancement property. Note that 

((BI + Ni)- 1 + d) ~*= ((Bi + Ni) _1 (I + (B! + Ni) Ox)) ~* (20) 

( = } (i + N^r 1 (b; + no - Bi + Bt 

= (I + N 1 1 y 1 ((Bi + - (I + N1O1) Bi) + Bi 

= (I + NiOi) -1 Ni + BJ 

= (Ni (Ni 1 + 1 )) _1 N 1 +Bi 

= (Ni 1 + 0!) _1 +Bi 

= Bi + n' 15 

where (a) and (6) follows from the fact that B^Oi = 0. Therefore, according to ( fTTI i the following property holds for the 
enhanced channel. 

(Bi + n^- 1 + ( M - i)(b; + N3)- 1 = M (B£ + N 2 ) -1 . 

Since ^ N 2 ^ N 3 then, there exists a matrix A such that N 2 = (I — A)N' 1 + AN 3 where A = (N 2 - N' 1 )(N 3 — N 1 ) _1 . 
Therefore, the above equation can be written as. 

(Bi+N' 1 )- 1 + ( M -l)(Bi+N 3 )- 1 = 



(I-A)(Bi + N 1 )+A(Bi+N s 



-1 



Let (I — A)(BJ + N x ) = aA(Bi + N 3 ) then after some manipulations, the above equation becomes 

-I + (^-l--)A = -^-I. (21) 
a a a + 1 
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The above equation is satisfied by a — —^j which completes the proportionality property. We can now prove the rate 



conservation property. The expression J — r^-i — 1 can be written as follow. 



l N 'i| 



Ni(B;+Ni) 



(22) 



HI 



(BJ+Ni-Bj) (BJ + Ni) 
HI 



I - B* (B* + N'J 

III 



|I-B*((B*+N x )-i 

(a) |I| 



Oi) 



I B* (BJ+Ni) 



|Ni| ' 

where (a) once again follows from the fact that BJOi = 0. To complete the proof of rate conservation, consider the following 
equalities. 

(b; 



B* + B* 



B'\ 



■n; 



(23) 



B5 (BI + N a ) 



w |b; 



N, 



|B*>N 2 | ' 

where (a) follows from the fact BiJjOa = 0. Therefore, according to (1221 . d23l . and the fact that N 3 = N3, the rate 
preservation property holds for the enhanced channel. To prove the optimality preservation, we need to show that (Bi,B 2 ) 
are also realizing matrices of an optimal Gaussian rate vector in the enhanced channel. For that purpose, we show that the 
necessary KKT conditions for the enhanced channel coincides with the KKT conditions of the original channel. The expression 
/i(Bj + B 2 + N 2 ) _1 can be written as follows 



(b; 



B? 



N. 



Mb; + b5 



N; 



io 2 



(24) 



H (b\ + B* 2 (l + B* 1 ^Na" 1 + io 




= M lB^+B* ^1 

( = } ^ (b; + b 2 (1 

= /x(b; + b^ (i + b;- 1 ^; 
= m(bj + b; + n 2 )- 1 

where (a) follows from the definition of N 2 and (6) follows from the fact that B 2 2 = 0. Therefore, according to d20l i, and 
the above equation, the KKT conditions of ( fTTb and ( f]~8T > for the original channel can be written as follows for the enhanced 
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channel. 



(B; + N;)- 1 + ( M - 1)(BJ + N3)- 1 = p(B; + N 2 )-\ (25) 
/i(B* + B* 2 + N 2 )-! = M (B* + B 2 + N3)- 1 + 3 2 . (26) 

where 3 2 h 0. Therefore, Jif (Bi, 2 , ^1,2,3) + ^R§{Bi }2 , N' 1)2)3 ) is maximized when B k = B£ for k = 1, 2. ■ 



We can now use Theorem[2]to prove that 1Z G (S, Ni^^) is the capacity region of the SADBC. We follow Bergman's approach 
[21] to prove a contradiction. Note that since the original channel is not proportional, we cannot apply Bergman's proof on 
the original channel directly. Here we apply his proof on the enhanced channel instead. 

Theorem 3: Consider a SADBC with positive definite noise covariance matrices (Ni, N 2 , N3). Let C(S,Ni. 2 ,3) denote 
the capacity region of the SADBC under a covariance matrix constraint S >- .Then, C(S, N1.2.3) = 1Z G (S, Ni i2 ,3)- 

Proof: The achievability scheme is secret superposition coding with Gaussian codebook. For the converse proof, we use 
a contradiction argument and assume that there exists an achievable rate vector R = (i?i,i?2) which is not in the Gaussian 
region. We can apply the steps of Bergman's proof of [21] on the enhanced channel to show that this assumption is impossible. 
Since R ^ 1Z G (S, Ni^.3), there exist realizing matrices of an optimal Gaussian rate vector BJ, B 2 such that 



for some b > 0. We know by Theorem that for every set of realizing matrices of an optimal Gaussian rate vector B^,B 2 , 
there exists an enhanced SADBC with noise covariance matrices N^N^ such that the proportionality and rate preservation 
properties hold. According to the rate preservation property, we have i?^(B^ 2 , Ni. 2 ) = i?^(B^ 2 , Nj 2 ), k = 1, 2. Therefore, 
the preceding expression can be rewritten as follows: 



Ri> ^f(B^ 2 ,Ni )2j3 ), 
R 2 >R$(Bl 2 ,N 1>2<3 ) + b, 



(27) 



R 1 > iZ?(B; )2 ,N 1)2(3 ) - #f(B£ )2 ,N' 1>2i3 ), 

R2> ^(B;, a ,N liai3 ) + b = i? 2 G (B^ 2 ,N' l!2i3 ) + b, 



(28) 



According to the Theorem [Tj R\ and R2 are bounded as follows: 



Ri< Myil u ) - h(z\u) - (ft(yi|x,u) - /i(z|x, u)) 



i? 2 < h(y 2 ) - h(z) - (h(y 2 \u) - h(z\u)) 



Let y ± and y 2 denote the enhanced channel outputs of each of the receiving users. As u — > y k — > y k forms a Markov chain 
for k = 1, 2 and z = z, then we can use the data processing inequality to rewrite the above region as follows: 




(29) 



h( yi \u) - h(z\u) - - (log IN'J - log |N 3 |)) 




(30) 



Now, the inequalities of d28T > and j29l have shifted to the enhanced channel. 
Since R\ > R G (Bi.^Nj 2 3 ), the inequality ( |29l means that 
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By the definition of matrix A and since y ± — > y 2 — > z forms a Morkov chain, the received signals z and y 2 can be written 
as z = y 1 + n and y 2 = y 1 + A^n where n is an independent Gaussian noise with covariance matrix N = N 3 — 
According to Costa's Entropy Power Inequality and the previous inequality, we have 



h(y 2 \u)-h(z \u) 




where (a) is due to the proportionality property. Using <f30b and the fact that R2 > i?2 ? (Bi,2, N x 2 3), the inequality (l30b 



means that 



h{y 2 )-h(z) > R 2 + h{y 2 \u) - h(z \u) > 

\ io g (B; + b; + n 2 ) - \ io g (B; + b* + n 3 ) 



On the other hand, Gaussian distribution maximizes h(x + n.2) — h(x + ris) (See [22]) and (BJ,B 2 ) satisfying the KKT 



In this section, we characterize the secrecy capacity region of the aligned (but not necessarily degraded) MIMO broadcast 
channel. Note that since the SAMBC is not degraded, there is no single-letter formula for its capacity region. In addition, the 
secret superposition of Gaussian codes along with successive decoding cannot work when the channel is not degraded. In [6], 
we presented an achievable rate region for the general secure Broadcast channel. Our achievable coding scheme is based on a 
combination of the random binning and the Gelfand-Pinsker binning schemes. We first review this scheme and then based on 
this result, we develop an achievable secret coding scheme for the SAMBC. After that, based on the Theorem [2] we provide 
a full characterization of the capacity region of SAMBC. 

A. Secret Dirty-Paper Coding Scheme and Achievability Proof 

In [6], we established an achievable rate region for the general secure broadcast channel. This scheme enables both joint 
encoding at the transmitter by using Gelfand-Pinsker binning and preserving confidentiality by using random binning. The 
following theorem summarizes the encoding strategy. The confidentiality proof is given in Appendix II for completeness. 

Theorem 4: : Let V\ and V% be auxiliary random variables and 57 be the class of joint probability densities P{v\, v^x, yi, 2/2, z) 
that factors as P(v\, V2)P(x\vi, V2)P(yi,U2, z\x). Let 1Zi(tt) denote the union of all non-negative rate pairs (Ri, R 2 ) satisfying 



for a given joint probability density ir 6 O. For the general broadcast channel with confidential messages, the following region 
is achievable. 



conditions of d26b . Therefore, the above inequality is a contradiction. 



IV. The Capacity Region of the SAMBC 



fli</(Vi;Yi)-I(7i;.Z), 



R2<I(V2;Y 2 )-I(V 2 ;Z), 



Ri + R2< HVliYl) + I{V 2 ;Y 2 ) - I(V U V 2 ; Z) - I{V X ; V 2 ), 




(32) 
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(Vf , V?) e 4"> 




Fig. 2. The Stochastic Encoder 

where conv is the convex closure operator. 

Remark 2: If we remove the secrecy constraints by removing the eavesdropper, then the above rate region becomes Marton's 
achievable region for the general broadcast channel. 

Proof: 1) Codebook Generation: The structure of the encoder is depicted in Fig|2] Fix P(vi), P(v 2 ) and P(x\vi, v 2 ). The 
stochastic encoder generates 2 n ( I ( Vl ' Yl '~ e > independent and identically distributed sequences t>" according to the distribution 
P( v i) — n™=i P{ v i,i)- Next, randomly distribute these sequences into 2 nRl bins such that each bin contains 2 n ( I ( Vi:Z )~ e ' ) 
codewords. Similarly, it generates 2™( / ( V2;Y2 ) _c ) independent and identically distributed sequences v 2 according to the distri- 
bution P(v 2 ) = n"=i P( v 2.i)- Randomly distribute these sequences into 2 nR2 bins such that each bin contains 2™( / ( v ' 2; ' z )~ e ) 
codewords. Index each of the above bins by i«i G {1, 2, 2 nRl } and w 2 6 {1, 2, 2 nR2 } respectively. 

2) Encoding: To send messages wi and w 2 , the transmitter looks for in bin w\ of the first bin set and looks for v% in bin 
w 2 of the second bin set, such that (w™, v 2 ) € Ai n \Py 1 y 2 ) where Ae (Pvi,v 2 ) denotes the set of jointly typical sequences t>™ 
and v 2 with respect to P(v\, v 2 ). The rates are such that there exist more than one joint typical pair. The transmitter randomly 
chooses one of them and then generates x n according to P(x n \v™, v 2 ) = Yi7=i P( x i\ v i,ii v 2,i)- This scheme is equivalent to 
the scenario in which each bin is divided into subbins and the transmitter randomly chooses one of the subbins of bin wi and 
one of the subbins of bin It then looks for a joint typical sequence (w™, v% ) in the corresponding subbins and generates 
x n . 

3) Decoding: The received signals at the legitimate receivers, an d y 2 , are the outputs of the channels P(y™\x n ) = 
n™=i P(lJi-i\ x i) and P(y 2 \x n ) = P{y2,i\%i), respectively. The first receiver looks for the unique sequence such that 
(vi,y 7 i) is jointly typical and declares the index of the bin containing as the message received. The second receiver uses 
the same method to extract the message w 2 - 

4) Error Probability Analysis: Since the region of (0 is a subset of Marton region, then the error probability analysis is the 
same as [3]. 

5) Equivocation Calculation: Please see Appendix A. ■ 
The achievability scheme in Theorem [4] introduces random binning. However, when we want to construct the rate region of 
f32) . it is not clear how to choose the auxiliary random variables V\ and V 2 . Here, we employ the Dirty-Paper Coding (DPC) 
technique to develop the secret DPC (SDPC) achievable rate region for the SAMBC. We consider a secret dirty-paper encoder 
with Gaussian codebooks as follows. 
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First, we separate the channel input x into two random vectors h>i and b2 such that 

t>i + b2 = x 

Here, t>i and b2 and vi and V2 are chosen as follows: 

bi~jV(0,B 1 ) ) 
b 2 ~AA(0,B 2 ), 
v 2 = b 2 , 
vi = bi + Cb 2 . 



(33) 



(34) 



where Bi = /?[bibj] y and B2 = /?[b2bj] y are covariance matrices such that Bi + B2 di S, and the matrix C is 
given as follows: 



B 1 (N, 



Bi 



(35) 



By substituting d34l into the Theorem [4] we obtain the following SDPC rate region for the SAMBC. 

Lemma 1: (SDPC Rate Region): Let S be a positive semi-definite matrix. Then the following SDPC rate region of an 
SAMBC with a covariance matrix constraint S is achievable. 



H SDPC (S,N lt2 , 3 ) = conv\ |J K SDPC (ir, S,N 1)2)3 ) 



(36) 



where J\ is the collection of all possible permutations of the ordered set {1,2}, conv is the convex closure operator and 
-]Z SDPC (TT,S,N lt 2,3) is given as follows: 



K 



SDPC 



(7r,S,N 1>2)3 ) = 



(R U R 2 ) \R k = RSPffifa Bi, 2 , Ni, a ,a) k = 1, 2 

s.t s - (Bi + b 2 ) y 0, Bi y 0, b 2 y o 



where 



^^( 7 r > B 1>2) N 1 , 2 , 3 ) = i 



log 



E^^B^+N, 
Note that for the identity permutation, 717, where TTj(k) = k we have 



\ log 



Er:i 1(fe) B. (i) +N 3 



7T- 1 (fc)-l 



Ei=l 



B 7r(i) + N 3 



K 



SDPC, 



7r/,S,N li2)3 ) =^ G (S,N li2 , 3 ) 



Proof: We prove the lemma for the case of identity permutation 717 = {1, 2}. This proof can similarly be used for the 
case that it = {2, 1}. According to the Theorem |4] we have, 

R t < min {I{Vv, Y x ) - /(VI; Z), I{V i; Yi) + I(V 2 ; Z) - I{V U V 2 ; Z) - I(V i; V 2 )} , 

< min {7(V1; F x ) - /(F x ; Z), J(V X ; Yi) - J(V X ; Z|V 2 ) - 7(V1; V 2 )} , 

< I(V 1 :Y 1 )-I(V 1 ;Z\V 2 )~I(V 1 ;V 2 ), 

R 2 <I{V 2 -Y 2 )-I{V 2 ;Z), (37) 

where (a) follows from the fact that /(Vi, V 2 ; Z) = /(V2; Z) + /(Vi; ZIV2) and (6) follows from the fact that /(Vi; Z\V 2 ) + 
/(li; 1^) = /(Z, l^; Vi) > /(Z; Vi). To calculate the upper-bound of R\, we need to review the following lemma which has 
been noted by several authors [23]. 
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Lemma 2: Let yi = t>i + b2 + ni, where bi, b2 and ni are Gaussian random vectors with covariance matrices Bi, 
B2 and Ni respectively. Let bi, b2 and ni be independent, and let vi = bi + Cb2, where C is an t x t matrix. Then 
an optimal matrix C which maximizes J(vi;yi) - /(vi;b 2 ) is C = Bi (Ni + Bi) . Further, the maximum value of 
^(vi;yi - 7(vi;b a ) is /(vi;yi|b a ). 

Now, using the above Lemma and substituting ( |34] i into ( 137) , we obtain the following achievable rate region when tt = 717. 



Ri< 

R2< 



log |Nj x (Bi + Ni)| - \ log |N 3 \Bx + N 8 ) 

|Bi +B 2 +N 2 | _ 1 |B! +B2 + N3I 
° g |B 1+ N 2 | 2 ° g |B 1+ N 3 | 



i + 



B. SAMBC- Converse Proof 

For the converse part, note that not all points on the boundary of 1Z SDPC (S, Ni.2,3) can be directly obtained using a single 
SDPC scheme. Instead, we must use time-sharing between points corresponding to different permutations. Therefore, unlike 
the SADBC case, we cannot use a similar notion to the optimal Gaussian rate vectors, as not all the boundary points can 
immediately characterized as a solution of an optimization problem. Instead, as the SDPC region is convex by definition, we 
use the notion of supporting hyperplanes of [13] to define this region. 

In this section, we first define the supporting hyperplane of a closed and bounded set. Then, we present the relation between 
the ideas of a supporting hyperplane and the enhanced channel in Theorem [5] This theorem is an extension of Theorem |2] to 
the SAMBC case. Finally, we use Theorem[5]to prove that 1Z SDPC (S, N1.2.3) is indeed the capacity region of the SAMBC. 

Definition 6: The set {R = (Ri, i?2)|7i-Ri + 72-R2 = b}, for fixed and given scalars 71, 72 and, b, is a supporting hyperplane 
of a closed and bounded set X C R"\ if 71-R1 + 72-R2 < b V(i?i,i?2) G X, with equality for at least one rate vector 
(Ri, R2) £ X. 

Note that as X is closed and bounded, max( fll ,R 2 )eX I1R1 + 72-R2, exists for any 71,72- Thus, we always can find a 
supporting hyperplane for the set X. As 1Z SDPC (S, N1.2.3) is a closed and convex set, for each rate pair of R = {RI, R2) ^ 
■fcSDPC '(SjNi^a) which lies outside the set, there exists a separating hyperplane {(Ri, R2)\jiRi + I2R2 = b} where 
7i > 0,71 > 0,6 > and 

71-Rx + 72 i?2< b V(Ri,R 2 ) G TZ SDPC {S, N w ) 
jiRl + j 2 R%> b 

The following theorem illustrates the relation between the ideas of enhanced channel and a supporting hyperplane. 

Theorem 5: Consider a SAMBC with noise covariance matrices (Ni,N a ,Na) and an average transmit covariance matrix 
constraint S y 0. Assume that {(i?i, R2)\jiRi+"f2R2 — b} is a supporting hyperplane of the rate region 1Z SDPC (717 , S, Ni ; 2,3) 
such that < 71 < 72, 72 > and b > 0. Then, there exists an enhanced SADBC with noise covariance matrices (N l5 N 2 , N 3 ) 
such that the following properties hold. 

1) Enhancement: 

N'^Ni, N 2 ^N 2 , N3-N3, N;^N 2 , 

2) Supporting hyperplane preservation: 

{(Ri, i?2)|7i-Ri + 72-R2 = b} is also a supporting hyperplane of the rate region 1Z G (S, N x 2 3 ) 
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Proof: To prove this theorem, we can follow the steps of the proof of Theorem [2] Assume that the hyperplane 
{(i?i,i?2)|7i-Ri + I2R2 = b} touches the region TZ SDPC (n I , S, Ni )2>3 ) at the pint (R^R^). Let Bi,B 2 be two positive 
semi-definite matrices such that B£ + B 2 S and such that 

Jef 2,pc (7ri,B; ia ,N 1 , a , 8 ) = R* k , k = 1,2 

By definition of the supporting hyperplane, the scalar b and the matrices (B£, B 2 ) are the solution of the following optimization 
problem: 

max 7l J2f DPC fa, B 12 , Ni.a.s) + j 2 RI DPC (tti, B 12) N 1i2 , 3 ) 

Bi .B2 

s.t Bi+B 2 ^S B k ^0 k = l,2 

We define the noise covariance matrices of the enhanced SADBC as ( fT9l ). Since for the permutation tt = 717 we have 
TZ sdpc (tti, S,N 1)2i3 ) = TZ G (S, Ni i2)3 ), the supporting hyperplane {(i?i, i? 2 )|7i-Ri + I2R2 = b} is also a supporting 
hyperplane of the rate region 1Z G (S, Nj 2 3)- ■ 
We can now use Theorem [5] and the capacity result of the SADBC to prove that 1Z SDPC (S, Ni j2 ,3) is indeed the capacity 
region of the SAMBC. The following theorem formally states the main result of this section. 

Theorem 6: Consider a SAMBC with positive definite noise covariance matrices (Ni, N 2 , N3). Let C(S,Ni, 2 ,3) denote 
the capacity region of the SAMBC under a covariance matrix constraint S >- .Then, C(S,Ni. 2 .3) = 1Z SDPC (S, Ni j2j 3). 
Proof: To proof this theorem, we use Theorem [5] to show that for every rate vector R° , which lies outside the region 
we can find an enhanced SADBC, whose capacity region does not contain R . As the capacity region of 
the enhanced channel outer bounds that of the original channel, therefore, R cannot be an achievable rate vector. 

Let R° = b e a rate vector which lies outside the region 1Z SDPC (S, Ni i2i3 ). There exists a supporting and 

separating hyperplane {(R\, i?2)|7i-Ri + 72-R2 = b} where 71 > 0, 72 > 0, and at least one of the 7fc's is positive. Without 
loose of generality, we assume that 72 > 71. If that is not the case, we can always reorder the indices of the users such that 
this assumption will hold. By definition of the region 1Z SDPC (S, N1.2.3), we have, 

n SDPC (^,s,N l!2 , 3 ) C ^ DPC (S,N l!2 , 3 ). 

Note that, as {(-Ri, -R 2 )|7i-Ri + 72^2 = b} is a supporting hyperplane of 1Z SDPC (S, Ni i2j3 ), we can wire, 

b= max 71-R1+72R2 

< max J1R1 + I2R2 = b. 

Furthermore, we can also write, 

71^+72^ >b>b. 

Therefore, the hyperplane of {(i?i, -R 2 )|7i-Ri + 72^2 = b } is a supporting and separating hyperplane for the rate re- 
gion TZ (717, S,Ni l2l 3). By Theorem [5] we know that there exists an enhanced SADBC whose Gaussian rate region 
7?. G (S, 2 3) lies under the supporting hyperplane and hence (_R°,i?2) ^ 7^ G (S,N 12 3)- Therefore, {RI^R^) must lies 
outside the capacity region of the enhanced SADBC. To complete the proof, note that the capacity region of the en- 
hanced SADBC contains that of the original channel and therefore, (i^iif) must lies outside the capacity region of the 
original SAMBC. As this statement is true for all rate vectors which lie outside 1Z SDPC (S, Ni^^), therefore we have 
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C(S,Ni, a ,3) C ^ SDPC (S,N l!2 , 3 ). However, TZ SDPC {S, Ni, 2 , 3 ) is the set of achievable rates and therefore, C(S, N li2i3 ) = 
R SDPC (S,N 1AS ). ■ 
With the same discussion of [13], the result of SAMBC can extend to the SGMBC and may be omitted here. The results of 
the secrecy capacity region for two receiver can be extended for m receivers as follows. 

Corollary 1: Consider a SGMBC with m receivers and one external eavesdropper. Let S be a positive semi-definite matrix. 
Then the SDPC rate region of ft SDPC (S,Ni,..., m ,Hi,..., m ), which is defined by the following convex closure is indeed the 
secrecy capacity region of the SGMBC under a covariance constraint S. 

K SDPC (S, Ni m , Hi m ) = conv J (J K sdpc (tt, S, N, ,„. H 1; ..., m ) I (38) 

Uen J 

where J| is the collection of all possible permutations of the ordered set {1, ...,m}, conv is the convex closure operator and 
H SDPC (ir,S,'Ni t ... tm ,'H. li ... im ) is given as follows: 

,„sz)/>r, c tv T , , I ( R i, R 2) |-Rfe = ^-im(' r . B i,.,m.Ni,.., I „,Hi,..., m ) k = l,...,m 



where 



s.t s - YZLi B * ^ °. B ^ °' * = m 

i?f ^ (tt, B, ,„. N 1; ..., m , H 1; ..., m )= i [log ■ 



H k 


(Er:; (fc) B. (i) ) 


H[ + N k 


H k ( 
H 3 ( 


Er:; (fe) B. (i) j 


(i) 


)Ht+I 
4 + N 3 


r 



(39) 



2' 
4 log 

V. Multiple-Input Single-Outputs Multiple Eavesdropper (MISOME) Channel 

In this section we investigate practical characterizations for the specific scenario in which the transmitter and the eavesdropper 
have multiple antennas, while both intended receivers have a single antenna. We refer to this configuration as the MISOME 
case. The significance of this model is when a base station wishes to broadcast secure information for small mobile units. In 
this scenario small mobile units have single antenna while the base station and the eavesdropper can afford multiple antennas. 
We can rewrite the signals received by the destination and the eavesdropper for the MISOME channel as follows. 

?/i = h\x + m, 

?/2 = hlx + n 2 , (40) 
z = H 3 x + n 3 , 

where hi and h 2 are fixed, real gain matrices which model the channel gains between the transmitter and the legitimate 
receivers. These are matrices of size t x 1. The channel state information again is assumed to be known perfectly at the 
transmitter and at all receivers. Here, the superscript f denotes the Hermitian transpose of a vector. Without lost of generality, 
we assume that m and n 2 are i.i.d real Gaussian random variables with zero means unit covariances, i.e., ni,n 2 ~ Af(Q, 1). 
Furthermore, we assume that n 3 is a Gaussian random vector with zero mean and covariance matrix I. In this section, we 
assume that the input x satisfies a total power constraint of P, i.e., 

Tr{E(xx T )} < P 

Before we state our results for the MISOME channel, we need to review some properties of generalized eigenvalues and 
eigenvectors. For more details of this topic, see, e.g., [24]. 
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Definition 7: (Generalized eigenvalue-eigenvector) Let A be a Hermitian matrix and B be a positive definite matrix. Then, 
(A, i/j) is a generalized eigenvalue-eigenvector pair if it satisfy the following equation. 

Atp = XBip 

Note that as B is invertible, the generalized eigenvalues and eigenvectors of the pair (A, B) are the regular eigenvalues 
and eigenvectors of the matrix B X A. The following Lemma, describes the variational characterization of the generalized 
eigenvalue-eigenvector pair. 

Lemma 3: (Variational Characterization) Let r(i/>) be the Rayleigh quotient defined as the following. 

,,s = ^Atp 

Then, the generalized eigenvectors of (A, B) are the stationary point solution of the Rayleigh quotient r(if>). Specifically, 
the largest generalized eigenvalue A max is the maximum of the Rayleigh quotient r(if>) and the optimum is attained by the 
eigenvector i/> max which is corresponded to A max , i.e., 

maxr(V) = t^sM^ = Amax 
V> V T m ax B Vmax 

Now consider the MISOME channel of (1401 . Assume that < a < 1 and P are fixed. Let define the following matrices for 
this channel. 

Ai,i= I + aPhjhJ, 

Bi, 1= i + apnln 3 

Suppose that (A(x,i) max , V'lmax) ^ s ^ e l ar g est generalized eigenvalue and the corresponding eigenvector pair of the pencil 
(Ai i,Bi i). We furthermore define the following matrices for the MISOME channel. 

A 2 2 = I + r 2 h 2 h 2 , 

l + oP|h 2 t/> lmax | 2 



B 



2,2 — 



I + (1 - a)PHl (i + «FH 3 ^ max </4 max H 3 ) 1 H 3 

Assume that (A( 2 ,2) maX ; V^max) i s tne largest generalized eigenvalue and the corresponding eigenvector pair of the pencil 
(A 2j2 ,B 2 . 2 ). Moreover, consider the following matrices for this channel. 

A 2!l =I+(l-a)Ph 2 h 2! 
B 2)1 =I+(l-a)PH|H 3 , 



aP 

Ai 2 = H 1 hih{. 

l + (l-a)P|hiV;3 max | 2 

B 12 = I + aPUl (i + (1 - a)PH 3 t/> 3max ^ max H 3 ) * H 3 , 

where we assume that (A( 2j i) max , V^rnax)' an( ^ CVi,2) m ax , / 04 max ) aie me largest generalized eigenvalue and the corresponding 
eigenvector pair of the pencils (A 2 i,B 2 .i), and (Ai 2 ,Bi. 2 ) respectively. The following theorem then characterizes the 
capacity region of the MISOME channel under a total power constraint P based on the above parameters. 

Theorem 7: Let Q MI SOME denote secl -g C y capacity region of the the MISOME channel under an average total power 
constraint P. Let Y[ be the collection of all possible permutations of the ordered set {1,2} and conv be the convex closure 
operator, then c MISOME is given as follows. 



c MISOME = J II U M ISOMER 
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where Ji MISOME ^ i s given as follows. 



n 



MISOME 



w= U n 



MISOME 



0<a<l 



where 1Z MISOME (n, a) is the set of all (Pi,P 2 ) satisfying the following condition. 



1 



Pfc< 75 [l°g^(fc,7T- 1 (fc))max] 



fc = 1,2. 

Proof: This theorem is a special case of Theorem [6] and corollary Q] First assume that the permutation ir = ttj = {1, 2}. 
In the SDPC achievable rate region of d39l ), we choose the covariance matrices Bi and B2 are as follows. 

B 1= aPVlma X </4ma X > 

B 2 = (l-«)PV 2max ^ max - 
In other words, the channel input x is separated into two vectors b x and b 2 such that 

x= bi + b 2 , 

bi= UlV'lroax! 
b 2 = U 2 ^2max- 

where u x ~ A/"(0, aP), u 2 ~ A/"(0, (1 - a)P), and < a < 1. Using these parameters, the region of TZ sdpc {tt i , S, N 1i2 , 3 ) 
becomes as follows. 



_ 1 
~~ 2 



log 
log 

log 



1 



log 

2 8 



1 + htBihi 

(l + aPh^ lmax Vi max h i) 

1 111 



I + HsBxHt 



-5 log 



(°) 1 n \ 1 + 

= 7j L l0 S A (ia)maxJ 



(41) 



where (a) is due to the fact that |I + AB| = |I + BA| and the fact that i/'J mas i/'imax = 1 Similarly for the R2 we have, 



R 2 < 



log 
log 

log- 



1 + 


h 2 (B 1 +B 2 


)h 2 




l + h 2 B!h 2 




1 + 


h 2 B 2 h 2 


1 


l + h 2 Bih 2 


~ 2 



•log 



I + H3 (B a +B 2 )H| 



I + HaBiHj, 



H 3 B 2 H 3 



</4 



I + H 3 BiH 3 

(l-a)Ph 2 h^ 



2 max ^ I+^pjhJ^ 



V>2 



2 max 



I + (1 - a)PH 3 I + aPHaVwM 



H T 

max 3 



H 3 W» ai 



= 2 [l°g A (2, 2) max] + 



(42) 
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Similarly, when tt — {2, 1}, in the SDPC region, we choose t>i = Uiip imax and b 2 = U2ip 3max - Then the SDPC region is 
given as follows. 



log 



log 



log 



1 + h! (Bi + B a ) hi 



1 + hlBahi 



I + H 3 (B!+B a )Hj 



I + H3B2H3 



1 + 



hlBihi 



l + h^Bahi 



log 



H3B1H3 



I + H 3 B 2 H 3 



^4 max [} + 1 + (1 _ Q)P | h t^ 3max | 2 J V 



4 max 



</>Lax (l + "PH 3 (i + (1 - a )PH 3 t/, 3max </, 3max H 3 ) ~' H 3 ) </> 



4 max 



— ^ [ lo § \l,2) max] + 



P 2 < 



and i? 2 is bounded as follows 

1 
2 
1 

2 



log 
log 



1 + h 2 B 2 h 2 



~2 log 



I + H 3 B 2 H 



log 



(1 + (1 - a)Ph 2 ^ 3max i/) 3max h 2 

4 m ax(l+(l-«)^h 2 h 2 )^ 



log 



(I + (l- a )PH 3 t/> 3max ^ 3max H 3 ) 



3 max 



^Lax(l+(1-«)^H 3 H 3 ) ^ 3max 



= ^ [ lo gA(2,l)max] + ■ 



Note that the eigenvalues A (;ife)max = A(; ife) max (a, P) and the eigenvector Vfcmax = ''/'fcmaxK f,fc = 1,2 are the 
functions of a and P. The following corollary characterizes the secrecy capacity region of the MISOME channel in high SNR 
regime. 

Corollary 2: In the high SNR regime, the secrecy capacity region of the MISOME channel is given as follows. 



lim c C MISOME = conv \ \J TZ^ 



MISOME 



(tt) 



where 



K 



MISOME 



(t = {1,2» 



(Pi,P 2 ),Pi < X - [logA max (hihJ.HjHs) 



K 



MISOME 



(7T={2,1}) = 



(Pl,P 2 ),Pl < 



log 



An^^xht.HjHs) 



+ 1 



log- 



(h 2 h 2 ,H 3 



3 H 3 



logA max (h 2 h 2 ,H 3 H £ 



where (A max (A ; , B), i/> imax ) denotes the largest eigenvalue and corresponding eigenvector of the pencil (A;,B) and b = 

l|H 3 t/' lmax || 2 ' HH 3 l/> 2max |p- 

Note that the above secrecy rate region is independent of a and therefore is a convex closure of two rectangular regions. 
Proof: We restrict our attention to the case that X(i_k) mM (a, P) > 1 for I, k = 1, 2 where the rates Pi and R 2 are nonzero. 
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First suppose that tt = 717 = {1,2}. We show that 

lim A (1 ,i) max (a,P)= A max (hih^H^Hg 



lim A( 2! 2)max(a,-P) = 



(h 2 h 2 ,H 3 H ; 



Note that since 



where 



for all P > we have, 



^(1,1) max 



(a,P) = 



> 1 



V'lmax^-P) = arg max 



l + aP\h[^ lmax (a,P)\ 
l + aP||H 3 Vi max («,P)|| 2 

l + o;P|hiVi(Q!,-P)| 2 



-[•0 1 :||V 1 ll 2 =i} l + a-P||H 3 Vi(a,-P)H 5 



1*4^1 max(«.^)| 2 > |H 3 ^ lmax («,P)|| : 



Therefore, An i i) max is an increasing function of P. Thus, 

A(i,i) ma x(o;,P)< 



|hl^ lmax (a,P)| 2 
j|H 3 t/* lmax (a,P)|| 2 



< A. 



(hihl,H 3 



3 H 3 



Since A max ( hih^, H 3 H 3 ) is independent of P we have 



lim A 



(1,1) max — A max 



Next, defining 



we have the following lower bound 



■0 1 (cxj) = arg max 



(hxh^H^Ha 



t„/. |2 



{V' 1 :||l/' 1 || 2 = l} ll H 3^l| 



js + a|h^ lmax (oo)| 5 



lim A (1)1)max (a,P)> lim -f 

p^oo i + a||H 3 Vi max (oo)|| 2 

= A max (h 1 h t 1 ,H 3 H 3 



As the lower bound and upper bound coincide then we obtain (03]). Similarly to obtain d44b note that since 

1 , (i-q)P|hty> 2max ( a ,p)| 2 



X ( n p\ _ l+aP|h^ lmax (q : P)P 

A( 2 , 2 )max(«,P) - - ■ (1 _ a) p||H,^ 9m . x(a ,P)|p > 1 



l+aP||H 3 ^ lmax (a,P)|| 



where 



for all P > we have, 



Therefore, we have 



(l-q)P|hty? 2max ( a ,P)| 2 
l+aP|h 2 l/' lmax (Q,P)| 2 



{^ 2 :||^ 2 || 2 = l} 1 , (l-a)P||H 3 ^,, 11 „ x ( a ,P)|| 2 
l+QP||H 3 l/' lmax (Q,P)|| 2 

(l- a )P|h^ 2max (q,P)| 2 ^ (l- a )P||H 3 t/, 2ma >,P)|| 2 
1 + aP|h 2 t/> lmax (a,P)| 2 1 + aP||H 3 ^ lmax (a,P)|| 2 

l h 2^ 2max (°°)| 2 

p lrm A (2 , 2)max (a,P)< ^ (oo)||a 



< 



H H 3^1 m ax(°°)l| 2 

A max (h 2 h 2 ,H 3 H 3 



(43) 
(44) 



(45) 
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where 

1*4^1 maxC 00 ) I 5 



b 



"02 (°°) = ar S max 



:l H 301 m a X (°o)l| 2 

l^2' i / , 2max| 2 

II 2 
max 1 1 



'{VV# 2 P=l} ll H 3</>2 

On the other hand we have the following lower bound 

I , (l-«)|h^ 2max (oo)| 2 

lim A( 2 ,2)max(«i-P) > 7. — t^tzt, r^r (46) 

•P^ 1 ^ ii (i-a)||H3'y 2max (oo)|| 2 

+ a!|H 3 ^ max (oo)P 

Note that < a < 1. It is easy to show that the right side of equation of ( l46b is a decreasing function of a and therefore the 
maximum value of this function is when a = 0. Thus we have, 

A max (h 2 h 2 ,H^H 3 ) 

lim A( 2 ,2)max(a,-P) > 7 

As the lower bound and upper bound coincide then we obtain d44b . When n = 2, 1, the proof is similar and may be omitted 
here. ■ 
Now consider the MISOME channel with m single antenna receivers and an external eavesdropper. Let x = Y^k=i where 
b k = Wfc^fcmax, u k ~ 7V(0, a k P), and ^feLi a fc = 1 - Assume that (Afc maXj V'femax) is the largest generalized eigenvalue and 
the corresponding eigenvector pair of the pencil 

(t, ^^.I + ^PHtfl + HaAHn^Hsl 



\ 1 + h^Ahk ' 

where A = (X^=i Q^(i)^VV(i) max^(j) max)- The following corollary then characterizes the capacity region of the 
MISOME channel with m receivers under a total power constraint P. 

Corollary 3: Let Y[ be the collection of all possible permutations of the ordered set {1, ...,m} and conv be the convex 
closure operator, then q mi Sume j s gj ven as follows. 

c MISOME = cony J (J K MMOMB (f) 

where 1Z MI {it) is given as follows. 

<r>MISOME r\ I I ^MISOME / „ „ \ 

0<Q fc <l.Er=i a*=l 

where 1Z MI (ir,ai, ...,a m ) is the set of all (i?i, ...,R m ) satisfying the following condition. 

•Rfc< ^ [log Afcmax] + , k=l,...,m. 

VI. Conclusion 

A scenario where a source node wishes to broadcast two confidential messages for two respective receivers via a Gaussian 
MIMO broadcast channel, while a wire-tapper also receives the transmitted signal via another MIMO channel is considered. 
We considered the secure vector Gaussian degraded broadcast channel and established its capacity region. Our achievability 
scheme was the secret superposition of Gaussian codes. Instead of solving a nonconvex problem, we used the notion of 
an enhanced channel to show that secret superposition of Gaussian codes is optimal. To characterize the secrecy capacity 
region of the vector Gaussian degraded broadcast channel, we only enhanced the channels for the legitimate receivers, and the 
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channel of the eavesdropper remained unchanged. Then we extended the result of the degraded case to non-degraded case. We 
showed that the secret superposition of Gaussian codes along with successive decoding cannot work when the channels are 
not degraded, we developed a Secret Dirty Paper Coding (SDPC) scheme and showed that SDPC is optimal for this channel. 
Finally, We investigated practical characterizations for the specific scenario in which the transmitter and the eavesdropper can 
afford multiple antennas, while both intended receivers have a single antenna. We characterized the secrecy capacity region in 
terms of generalized eigenvalues of the receivers' channels and the eavesdropper channel. In high SNR we showed that the 
capacity region is a convex closure of two rectangular regions. 



A. Equivocation Calculation 

The proof of secrecy requirement for each individual message (0 and ( fTOb is straightforward and may therefore be omitted. 
To prove the requirement of (fTTT i from H (Wi, W 2 \Z n ), we have 



nR el2 = H{W u W 2 \Z n ) 

= H{W!,W2,Z n )-H(Z n ) 

= H(Wi,W 2 , Vr, V 2 n , Z n ) - H{V?, V?\W U W 2 , Z n ) - H(Z n ) 

= H{W U W 2 , Vr,V 2 n ) + H{Z n \W u W 2 ,V?, V 2 n ) - H{V?, V 2 n \W u W 2 , Z n ) - H(Z n ) 

> H{W U W 2 , Vr,V 2 n ) + H{Z n \Wi, W 2 ,Vr, V 2 n ) - ne n - H(Z n ) 

( = } H(W U W 2 , V?, V 2 n ) + H{Z n \Vr, V 2 n ) - ne n - H{Z n ) 



> H(Vf, V 2 n ) + H{Z n \Vl\V 2 n ) - ne n - H(Z n ) 

= H(Vn + H(V 2 n ) - I(V?; V 2 ) - I(Vr,V 2 n ; Z n ) - ne n 

> I{V?; Y?) + I{V 2 n ;Y 2 n ) - I{V?; V 2 n ) - I{V?,V?; Z n ) - ne, 

(e) 

> nRi + nR 2 — ne n , 



where (a) follows from Fano's inequality, which states that for sufficiently large n, H (V™, V 2 n \Wx, W 2 , Z n ) < h{PwJ) 
+nP™ e R w < ne n . Here P™ e denotes the wiretapper's error probability of decoding (w™, v 2 ) in the case that the bin numbers 
wi and w 2 are known to the eavesdropper. Since the sum rate is small enough, then _P™ e — > for sufficiently large n. (b) follows 
from the following Markov chain: {W u W 2 ) -» {V?,V£) -> Z n . Hence, we have H(Z n \W u W 2 , V?, V 2 n ) = H(Z n \Vf, V 2 n ). 
(c) follows from the fact that H(Wi,W 2 ,V{ 1 ,V 2 n ) > H(V?,V 2 n ). (d) follows from that fact that H(V{ 1 ) > IiVfiY?) and 
H{V 2 n ) > I(V 2 n ;Y 2 n ). (e) follows from the following lemmas. 

Lemma 4: Assume V", V 2 and Z n are generated according to the achievablity scheme of Theorem |4] then we have, 



Proof: Let A^(Pv- L ,v 2 ,z) denote the set of typical sequences (Vj™, V 2 n , Z n ) with respect to Pv x ,V2,Z> and 
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I(V?, V 2 n ; Z n ) < nI(V 1} V 2 ; Z) + n5 v 



I{V?;V 2 n ) <n/(F i; y 2 ) + n5 2 , 




1, (V?,V?,Z n ) $ A?(Pv uVa , z ); 
0, otherwise, 
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be the corresponding indicator function. We expand I(V™, V 2 ; Z n ) as follow, 

I(Vr, V 2 n ;Z n )<I{V?,V?,Q;Z n ) (47) 

= I(Vl\V 2 n ;Z n X)+I((;Z n ) 
1 

= E p (c= mvr, v 2 n - z", c = j) + i(c, z n )- 

3=0 

According to the joint typicality property, we have 

P(( - 1)I(V?, V 2 n ; Z n \C = 1)< nP{(V?, V 2 n , Z n ) £ A?(P Vl ,v 2>z )) log \\Z\\ (48) 

<ne„log||Z||. 

Note that, 

7(C; Z n ) < H(C) < 1 (49) 
Now consider the term P(( = 0)I(V", V 2 n ; Z n \( = 0). Following the sequence joint typicality properties, we have 

P(C = 0)I(Vr, V 2 n ; Z n \C = 0)< I(Vr, V 2 n ; Z n \( = 0) (50) 

P(y?, v 2 n , z n ){ log p(yr, vi\ z n ) - log p(y?, v 2 n ) 

-logP(Z n )), 

< n [-H(V U V 2> Z) + H(V lt V 2 ) + H(Z) + 3e n ] , 
= n[I(V 1 ,V 2 ;Z)+3e n }. 

By substituting d48l . j49l , and (l50l l into d47] >, we get the desired reasult, 

I(Vr, V 2 n ; Z n )< nI{V x , V 2 ; Z)+n ^e„ log \\Z\\ + 3e„ + ^ , (51) 
= nI(V 1 ,V 2 ;Z) + nS ln , 

where, 

Sm = e«log \\Z\\ + 3e„ + -. 

n 

Following the same steps, one can prove that 

I(V?; V 2 n ) < nI(Vi;V 2 ) + nS 2n . (52) 

■ 

Using the same approach as in Lemma [4] we can prove the following lemmas. 

Lemma 5: Assume V™ , Y™ and Y 2 are generated according to the achievablity scheme of Theorem [4] then we have, 

I(V?;Y?) <nI(Vi;Yi)+nS 3n , 

I(V 2 n ;Y 2 n )<nI(Vr,Z)+nS An . 
Proof: The steps of the proof are very similar to the steps of proof of Lemma @] and may be omitted here. ■ 
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